Experimental Realization of Deutsch's Algorithm in a One-way Quantum Computer 



M. S. Tame 1 , R. Prevedel 2 , M. Paternostro 1 , P. B6hi 2 Q M. S. Kim 1 , and A. Zeilinger 2 ' 3 
1 School of Mathematics and Physics, Queen's University, Belfast BT7 INN, United Kingdom 
2 Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090 Vienna, Austria 
3 Institute for Quantum Optics and Quantum Information (IQOQI), 
Austrian Academy of Sciences, Boltzmanngasse 3, A-1090 Vienna, Austria 
(Dated: February 1, 2008) 

We report the first experimental demonstration of an all-optical one-way implementation of 
Deutsch's quantum algorithm on a four-qubit cluster state. All the possible configurations of a 
balanced or constant function acting on a two-qubit register are realized within the measurement- 
based model for quantum computation. The experimental results are in excellent agreement with 
the theoretical model, therefore demonstrating the successful performance of the algorithm. 
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The increasing interest in topics of quantum infor- 
mation processing (QIP) and quantum computation has 
stimulated considerable efforts in the realization of quan- 
tum hardware based on various experimental settings. 
These efforts have resulted in the realization of one and 
two-qubit logical gates [l| , even though the networking of 
these basic building blocks is still far from being practi- 
cal. Nevertheless, investigations in this direction, both at 
the experimental and theoretical level are vital for the ad- 
vancement of QIP. The ultimate aim is the realization of 
multi-qubit quantum algorithms able to outperform their 
classical analogues [HQ- I n this context, the implemen- 
tation of few-qubit quantum algorithms represents a step 
forward in the construction of working processors based 
on quantum technology d, 0] . 

Very recently, a radical change of perspective in the de- 
sign of quantum computational protocols has been pro- 
posed and formalized in the "one-way" model [B} . Here, 
computation is not performed by inducing a sequence of 
logical gates involving the elements of a quantum regis- 
ter, as in the quantum circuit model [lj. In the one-way 
case, a multipartite entangled state, the cluster state, is 
used as a resource for running a "program" represented 
by single-qubit measurements, performed in order to sim- 
ulate a given computational task 5] . This new paradigm 
for quantum computation, which limits the amount of 
control one needs over a register to the ability of perform- 
ing single-qubit measurements, has raised an enormous 
interest in the physical community. It has triggered inves- 
tigations directed toward a better understanding of the 
model 0] and also its practical applications 0, S| ■ The 
efforts produced so far have culminated in the experimen- 
tal demonstration of the basic features of the model, the 
realization of a two-qubit quantum search algorithm 7] 
and the theoretical proposal for a measurement-based re- 
alization of a quantum game [9j. The one-way model is 
also helping us to understand the paramount role of mea- 
surements in the quantum dynamics of a system. 

Here, we report the first experimental demonstration 
of a one-way based implementation of Deutsch's algo- 



rithm 3]. It represents a simple but yet interesting in- 
stance of the role that the inherent parallelism of quan- 
tum computation plays in the speed-up characterizing 
quantum versions of classical problems. We have used 
an all-optical setup, where the construction of cluster 
states has been successfully demonstrated 0, H| • Negli- 
gible decoherence affecting qubits embodied by photonic 
degrees of freedom ensure the performance of the proto- 
col in a virtually noise-free setting. Although Deutsch's 
algorithm has been implemented in a linear optical setup 
before [To} , our protocol represents its first realization in 
the context of one-way quantum computation. It is based 
on the use of an entangled resource locally equivalent to 
the cluster state used previously for performing a two- 
qubit search algorithm [tJ and reinforces the idea of the 
high flexibility of cluster resources. We show that four 
qubits in a linear cluster configuration are sufficient to 
realize all the possible functions acting on a logical two- 
qubit register. Two of these result from the application 
of an entangling gate to the elements of the register. In 
principle, this gate can be realized by inducing an inter- 
action between the photonic qubits. In our cluster state 
approach, the required entangling operations are real- 
ized by using the entanglement in the cluster resource 
and the nonlinearity induced by detection. There is no 
need for engineering it in a case by case basis (Tpj | , which 
is a very important advantage. The density matrix of 
the logical output qubits for the functions show excellent 
performance of the algorithm in our setup. 

Model- The generalized version of Deutsch's algo- 
rithm, also known as the Deutsch-Josza algorithm llj . 
takes an iV-bit binary input x and allows one to distin- 
guish two different types of function f(x) implemented 
by an oracle. A function is constant if it returns the 
same value (either or 1) for all possible inputs of x and 
balanced if it returns for half of the inputs and I for 
the other half. Classically one needs to query this oracle 
as many as 2 Ar ~ 1 + 1 times in some cases. However the 



quantum version requires only one query in all cases 11 1 . 
In the two-qubit version the algorithm implements 
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FIG. 1: Network diagrams for the black boxes in Deutsch's 
algorithm. We have BB(i)=fl ® 1, BB(ii)=l ® a x , 
BB(iii)=CNOT and BB(iv) = (l ® cr^CNOT (CNOT denotes 
a Control-NOT gate). 



the oracle as a function / on a single query bit x using 
an input ancilla bit y. The applied unitary operation 
is given by \x) \y) — > |x) |y © f(x)}. Preparing the input 
state as |+) |-), where |±) = (|0)±|l))/v / 2 and {|0> , |1>} 
is the single-qubit computational basis, the oracle maps 
the state to (l/V2)[(-l) /{0) |0) + (-l)/« |1)] |-). By 
measuring the query qubit in the {|±)} basis, one can 
determine which type of function f(x) corresponds to. 
If f(x) is balanced (constant), the query qubit is always 
|— ) (|+)). Thus, only one query of the oracle is necessary, 
compared to two in the classical version. 

The action of the above oracle is either preset or dic- 
tated by the outcome of another algorithm. In order to 
implement all possible configurations that it might take 
in the two-qubit version, we must be able to construct 
them using a combination of quantum gates. In Fig. [T] 
we show all possible oracles in terms of their quantum 
network. By describing each as a "black box" , one can 
see that all four black boxes (BB(i)-(iv)) implement their 
respective oracle operation. In order to carry out the 
algorithm using these quantum gates, we use a cluster 
state and carry out one-way quantum computation on it 
by performing a correct program of measurements. No 
adjustment to the experimental set-up is necessary. 

Given a cluster state, there are two types of single- 
qubit measurements that allow a one-way quantum com- 
puter to operate. First, by measuring a qubit j in the 
computational basis it can be disentangled and removed 
from the cluster, leaving a smaller cluster state of the re- 
maining qubits. Second, in order to perform QIP, qubits 
must be measured in the basis Bj(oi) = {|a+) ■ , ja—)^}, 

where \a±) j = (|0) ± e ta \l))j/\/2 {a e R). Choosing 
the measurement basis determines the rotation R z (a) — 
exp(— iaa z /2), followed by a Hadamard operation H = 
(a x +a z )/^/2 being simulated on an encoded logical qubit 
in the cluster residing on qubit j {o x ,y,z are the Pauli ma- 
trices). With a large enough cluster, any quantum logic 
operation can be performed with a proper choice for the 



Experimental implementation.- For the entangled re- 
source, in an ideal case, the following four-photon state 
is produced by means of the set-up shown in Fig. O (a) 
|$ c ) = (1/2)(|0000) + |0011) + 1 1100) - |Hll))i234 with 
(l-Oj) embodied by the horizontal (vertical) polar- 
ization state of one photon populating a spatial mode 
j = 1,..,4. The preparation of the resource relies on 
postselection: a four-photon coincidence event at the 
detectors facing each spatial mode witnesses the prepa- 
ration of the state. This state is locally equivalent to 
a four-qubit linear cluster state \$n n ) (the local oper- 
ation being H± ® I2 Cg> I3 ® H4). The experimentally 
produced state g is verified by means of a maximum- 
likelihood technique for tomographic reconstruction (l3j 
performed over a set of 1296 local measurements 0], each 
acquired within a time-window of 500 s. This provides 
information about the overall quality of the experimental 
state on which the algorithm is performed. We have used 
all the possible combinations of the elements of the mu- 
tually unbiased basis {|0) , |1) , |+> , |-) , \R) , \L)} 3 with 
|±) • embodied by the polarization state at ±45° and 
\L/R)j = (|0) ± i\l))j/y/2 corresponding to left and 
right-circularly polarized photons. This over-complete 
state tomography has the advantage of providing a more 
precise state estimation and significantly smaller error 
bars Q. The reconstructed density matrix of g is shown 
in Fig. [21(c) & (d) and has a fidelity with the ideal state 
of F = ($ c |(?|$ c ) = 0.62 ± 0.01. The error bar was 
estimated by performing a 100 run Monte Carlo simula- 
tion of the whole state tomography analysis, with Poisso- 
nian noise added to the count statistics in each run [13l | . 
Obtaining a higher fidelity is limited by phase instabil- 
ity during the lengthy process of state tomography and 
non-ideal optical elements. However, it is well-above the 
limit F = 0.5 for any biseparable four-qubit state 14 1 
and demonstrates the presence of genuine four particle 
entanglement. 

In order to perform Deutsch's algorithm on the cluster 
resource |$c)j we have used a specific set of measurement 
bases for the qubits in each black box case. In Table U 
we provide these basis sets (BB C ) and feed-forward (FF) 
operations used to carry out the black boxes on |$ c ) and 
also \$>un) (BB basis sets). As BB(ii) and BB(iv) are 
obtained from BB(i) and BB(iii) by using alternative FF 
operations (corresponding to adaptive measurements on 
the output qubits fl2|). in what follows we explicitly de- 
scribe BB(i) and BB(iii). Fig. [2] (b) shows the in-out 
logical states of the algorithm, where the logical input 
state corresponding to |x) = |+) is encoded on qubit 
1. The state \y) = \— ) will be encoded on qubit 3 by 
measuring qubit 4 in the B^n) basis during the imple- 
mentation of the algorithm (described next). This gives 
|a:)|y) = (l®fl i (7r))|+)|+). 

Qubit 2 in \$u n ) plays the pivotal role of the oracle as 
it performs a two-qubit gate on the logical input states 
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|x) and \y). For BB(i), measuring qubit 2 in the compu- 
tational basis disentangles it from the cluster and |$M n ) is 
transformed into |±) 1 (l/y/2)(\0) \+) ± |1) |-)) 34 (+ (-) 
for outcome |0) 2 (|1) 2 ))- The effective operation per- 
formed by this choice of the oracle's measurement basis 
is 1 <8> 1. By including the H operation applied to the 
input state \y) from the measurement of qubit 4, the over- 
all computation results in (1 ® ® HR z {-k)) |+) |+) 
which is equivalent to \x) \y® f(xj) = (U ® 1) |+) j— ) 
up to a local rotation H on physical qubit 3, applied 
at the FF stage. Qubits 1 and 3 can now be taken 
as the output |a;) \y © /(#))• For BB(iii), upon measur- 
ing qubit 2 in the B2(ir/2) basis, the oracle applies the 
gate 0R z (tt/2) <g> j R z (tt/2))CPHASE on |x) and \y) (see 




FIG. 2: (a): Experimental setup. An ultraviolet pump-laser 
performs two passages through a nonlinear Beta-Barium- 
Borate crystal (BBO) aligned to produce entangled photon 
pairs of the form (|00) - |ll}) ab /v > 2 and (|00) + |ll>)erf/V2- 
Compensators (Comp) are half-wave plates (HWP) and BBO 
crystals used in order to counteract walk-off effects at the 
BBO. By considering the possibility of obtaining a double- 
pair emission into the same pair of modes and the action of 
the polarizing-beam splitters (PBS's), the four terms enter- 
ing |$c) are obtained and their amplitudes and respective 
signs adjusted [3] with an additional HWP in mode a. The 
algorithm is executed by using quarter- wave plates (QWPs), 
HWP's, PBS's and photocounter pairs {Dj, Dy} for the per- 
formance of polarization measurements in arbitrary bases of 
the photons in mode j. (b): Sketch of the cluster-state con- 
figuration. Qubit 1 embodies the logical input for jx) and its 
output. Qubit 4 (3) is the logical input (output) for \y), which 
is always found to be |— ) 3 . (c) & (d): Real and Imaginary 
plots respectively of the reconstructed experimental density 
matrix g. 





Measurement basis 


BB(i) 


{B 1 (0),{|0) 2 ,|1) 2 },{|0)3,|1) 3 },B4W} 


BB c (i) 


{{10)! , \1) 1 },{\0) 2 , H>2>, {l°>3 ,|1> 8 M|1>4 . I°) 4 }} 


BB(iii) 


{S 1 ( 7 r/2),B 2 (7 r /2),{|0)3 J |l) 3 } ) B 4 (7 r )} 


BB c (iii) 


{Bi(3jt/2), B 2 (7t/2), {|0> a , |l),j, {|1) 4 , |0> 4 » 



TABLE I: Measurement bases for the black boxes. The FF 
operations are (a s x 2 ) 1 (a^) 3 for BB c (i) and (ap ® S4 )i(ctJ 4 ) 3 
for BB c (iii). Here, Sj is (1) if the outcome is (\ a -)j) 
on qubit j. 



Tame et al. in [6]), where CPHASE shifts the relative 
phase of the state |1) |1) by ir. This gives the com- 
putation \x)\y®f(x)) = CNOT|+)|-) = {R z {-k/2)® 
i? z (7r/2))CPHASE(l ® HR z {n)) \+) |+) up to local rota- 
tions R z (—it/2) <E> H R z (— 7t/2) on qubits 1 and 3, ap- 
plied at the FF stage. The measurements and outcomes 
of qubits 1, 3 and 4 constitute the algorithm. The ad- 
ditions to the FF stages described above, together with 
the measurement of qubit 2 should be viewed as being 
carried out entirely by the oracle. 

The results of our experiment are shown in Fig. [3l 
where we fully characterize the output states of our quan- 
tum computer by repeating the algorithm a large number 
of times. A single run of the algorithm (measuring the 
output qubit 1 in a specific basis only once) is sufficient 
in our setup to carry out the quantum computation with 
success rates as large as 90% (78%) for BB(i) (BB(iii)). 
However, repeating it several times allows us to verify 
the density matrix for the quantum state of qubits 1 
and 3 reconstructed through a maximum likelihood tech- 
nique [13|. Although only the logical state residing on 
qubit 1 provides the outcome of the algorithm, it is useful 
for the characterization of the quantum computer's per- 
formance to also determine the state residing on qubit 3. 
Ideally, the joint state of qubits 1 and 3 should be the 
product state \x)\y ® f{x)). By obtaining both correct 
logical output states, we can confirm that the algorithm 
will run correctly if included in a larger protocol. Fig. [3] 
shows the output density matrices for BB(i) and BB(iii). 
Both the no-feed-forward (no-FF) and FF situations are 
shown. In the latter case, the state of the output qubits is 
corrected from the randomness of the measurements per- 
formed on the physical qubits 2 and 4. From the previous 
analysis, we know that the expected outcome from a sin- 
gle run, when a constant (balanced) function is applied is 
| + , — ) 13 (| — , — ) 13 ). Evidently, the reconstructed density 
matrices, both in the FF and no-FF cases, show a very 
good performance of the algorithm when compared with 
the theoretical expectations. The real parts are domi- 
nated by the correct matrix elements and no significant 
imaginary parts are found. Quantitatively, the fidelity 
with the desired state in the case of a constant (balanced) 
function is found to be as large as 0.90±0.01 (0.78±0.01) 
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FIG. 3: The output density matrices for cluster qubits 1 and 3 
when BB(i) and BB(iii) are implemented. Panels (a) and (c) 
show the real parts of the two-qubit density matrix elements 
as obtained from a maximum likelihood reconstruction for 
the no-FF cases of BB(i) and BB(iii) respectively. Panels (b) 
and (d) show the corresponding plots for the FF case, due to 
the randomness of measurement outcomes for qubits 2 and 4. 
In all four cases the imaginary parts are zero in theory and 
negligible in the experiment (average values < 0.02). 

for the FF case and 0.82 ± 0.01 (0.63 ± 0.01) for the no- 
FF one. Moreover, no entanglement is found in any of 
the joint output states, as witnessed by the negativity of 
partial transposition criterion (lBj . The small admixture 
of the undesired |+, — ) 13 to the expected | — ,— ) 13 state 
when a balanced function is applied (Fig. [3] (c)) is due 
to the non-ideal fidelity of the experimental cluster state 
with This is more pronounced for BB(iii) than for 

BB(i), where the measurement basis of qubit 2 breaks 
the channel between |a;) and \y) resulting in a protocol- 
dependent noise-inheritance effect for imperfect cluster 
states (see Tame et al. in 0]). 

Remarks.- We have designed, demonstrated and char- 
acterized the performance of the first experimental re- 
alization of Dcutsch's algorithm on a four-qubit clus- 
ter state. Our experiment is one of the few quantum 
algorithms entirely implemented utilizing the one-way 
model 0, • The agreement between the experimental 
data and theory is excellent and only limited by the over- 
all quality of the entangled resource in the experiment. 
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